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Abstract 

We apply the displaced-oscillator variational ansatz to the Caldeira-Leggett model for a quantum particle in a one-dimensional 
box described by a tight-binding chain. We focus on the case of an Ohmic environment and study the phase diagram for 
different chain lengths. At zero temperature there is a phase transition to a localized phase when the number of sites is even. 
At finite temperature, a transition from a coherent to an incoherent regime is predicted for all the chain lengths considered. 
Finally, the results are compared to those obtained with numerical techniques. 
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1. Introduction 

Open quantum systems have been widely explored 
in the literature. From the very beginning of Quan- 
tum Mechanics it was realized that every real quan- 
tum system is, to some extent, coupled to an environ- 
ment. Hence, the investigation of the way this coupling 
can affect the properties of the system is crucial to un- 
derstand the experimental information. Moreover, the 
study of open quantum systems has proved very useful 
to analyze the quantum origin of dissipation [1]. The 
concept of environment has become central in the un- 
derstanding of how the classical world emerges from 
the quantum world [2] . 

In order to analyze open quantum systems, simple 
models of bath and system-bath coupling have been 
proposed. A bath composed of an infinite sum of har- 
monic oscillators is known to provide a sufficient de- 
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scription to understand a wide spectrum of properties 
of those systems [3] . As Caldeira and Leggett realized 
in their seminal paper on quantum dissipation [4], to 
reproduce the classical equations of motion of a particle 
with friction it is enough to consider a linear coupling 
to its position. Since then, the Caldeira-Leggett model 
has become a fruitful tool to understand basic proper- 
ties of the influence of the environment on a quantum 
particle in both dynamical and equilibrium contexts. 

In particular, the case of a one-dimensional quan- 
tum particle coupled to an environment exhibits a rich 
physical behavior. Three simple situations can be con- 
sidered. Firstly, that of a free particle coupled to the 
bath. In Ref. [5], it was shown that the effect of the 
bath is an effective renormalization of the mass of the 
particle, that can be understood physically as the par- 
ticle dragging a cloud of bath oscillators. This is actu- 
ally equivalent to the interpretation found in the po- 
laron problem [6] , where a related -but different- model 
is considered. The second case of interest is the parti- 
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cle in a periodic potential, studied in [7,8]. A remark- 
able difference with the first case is the existence of 
a quantum phase transition at zero temperature, in 
which the particle becomes localized. Finally, a parti- 
cle in a confined region of space can be considered (the 
so-called quantum particle in a box) . The lack of sym- 
metries of this model makes it difficult to apply most of 
the standard techniques available for other problems. 
Recently, a discretized version of the particle-in-a-box 
problem has been studied with the numerical renor- 
malization group (NRG) , performing an analysis of the 
zero-temperature phase diagram [9] . It was found that 
the coupling to the environment gives rise to a narrow- 
ing of the spatial density distribution of the particle, 
until a critical value of the coupling is reached where a 
localization phase transition occurs. 

In this paper we apply a variational ansatz to study 
the zero and finite temperature phase diagram of the 
dissipative quantum particle confined to a finite chain. 
The variational ansatz has been extensively used to 
analyze the celebrated spin-boson model (see [10] for 
a general description), which actually turns out to be 
the two-site limit of the model considered here. Here 
we will analyze the case of M sites, complementing and 
clarifying previous results at zero temperature, already 
reported in [9,11], and adding the finite temperature 
calculation. 

The paper is organized as follows: Section 2 describes 
the model, pointing out previous results obtained in the 
study of its equilibrium phase diagram at zero temper- 
ature. Section 3 contains the main results of the paper: 
the variational approach is introduced, and results are 
obtained for zero and finite temperature. Finally, sec- 
tion 4 contains a general discussion and a comparison 
with the NRG results. 



2. The model 

We analyze a finite tight-binding chain coupled to a 
dissipative bath in the way Caldeira and Legget envis- 
aged: 



T~t = Hkin + Hbath + Hint + Hct 

M 

Hkin = -t^2(\m)(m + l\ + \m + l)(m\) 

m 

Hbath = ^ kb{b k 

k<u c 

Hint = M J2 ^( b l + 6fc ) 

Hct = A 2 g 2 ]T 1 (1) 

k<u c 

The first term is the kinetic term of the tight-binding 
chain, m labelling the different M sites, and t being the 
hopping between sites. By taking t = h/2ma 2 , with m 
the mass of the particle and a the space between sites, 
in the limit a — > 0, M — > oo and L = Ma constant, we 
recover the continuum model of a confined quantum 
particle. The second term is that of a gapless bath of 
harmonic oscillators. The distribution of momentum of 
the oscillators can be taken continuous below a certain 
cutoff lu c . The third term is the linear coupling of the 
bath to the position q — ~}2 m (m — mo)\m){m\ of the 
particle, mo being a label for the center of the chain. 
The coupling is characterized by the bath spectral func- 
tion J(uj) = 7rA 2 X]fc<u> k5(u> ~ m our particular 
case, we will limit the discussion to an Ohmic bath, 
which fulfills the constraint J(u>) = 2ira\Lj\8(Lj c — lu), 
with a = \ 2 /4tt. Finally, the last term is a counter- 
term required to ensure a homogeneous dissipation 
along the chain despite the coupling to the position, 
which is sufficient in equilibrium contexts [12]. Actu- 
ally, that the coupling does not privilege any site can 
be seen by performing a unitary transformation on the 
Hamiltonian, U = e~ Xq '^k {h k~ bk ^ / ^ . As a result, 

M 

« = kblb k -tJ2[\^{m+l\e- X ^ ib ^ bk)/VK +H.c.] 

k<u> c m 

(2) 

Notice that this transformation diagonalizes exactly 
the Hamiltonian when t — 0. An important remark 
concerning this model is that it only preserves parity 
symmetry. 

As mentioned in the introduction, this model has 
been studied in Ref. [9] by using a version of the NRG 
technique conveniently suited to deal with bosonic 
baths [13]. Two remarkably results were found, (i) 
For small hopping (as compared to the cutoff, t <C 1), 
and couplings smaller than a certain critical coupling 
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of order a ~ 1, the coupling to the bath increasingly 
confines the particle to the center of the chain; i.e., 
despite of the homogeneous coupling, the boundary 
conditions are responsible for a non-homogenenous 
effect, (ii) For couplings larger than the critical value 
a new phase of the model was found, in which the 
particle gets localized at one single site, breaking the 
parity symmetry. However, not every site is permitted, 
only two sites localized symmetrically with respect to 
the center but never at the edges. In a naive analysis 
of the model an equal probability of getting localized 
at every site of the chain would have been expected. 
Thus, there is a second non-homogeneous effect com- 
ing from the presence of boundaries that shows up in 
the localized phase. 



3. Variational approach 

We will use the approach proposed by Silbey and 
Harris [14] to study the spin-boson model, and later ex- 
tended and discussed in a variety of papers [15,16,17]. 
The idea is to generalize the unitary transformation 
used in the previous section to get Hamiltonian (2) by 
the introduction of a variational parameter f k : 

U = e«E»<A/ fc ><*I- fc *> (3) 

As this kind of transformation displaces the oscillators 
from their unperturbed equilibrium positions, such a 
variational approach is conventionally referred to as the 
displaced oscillator ansatz. The transformed Hamilto- 
nian reads: 

H = UHU 1 * = Ho + H ICS 



Ho = -i rC n ^(|m)(m + 1| + H.c) + J2 kb l bk + 



9Q 



k<oj c 



Hrcs = V+ \ m + !)H + V -J2 l m X m V °1 

m m 

(4) 

where we have defined the potential strength g = 
J2 k (X + -^j) 2 and a renormalized hopping given 

by t rcn = (te~^k ifk/k)(b l- bk) ) T , with the ther- 
mal average taken with respect to Ho- The lat- 
ter can be worked out analytically, yielding t Icn = 
t exp[— i J2k W c °th(^r)]- Such a separation of the 



Hamiltonian ensures that the operators contained in 
the residual term H IC s, 



V+ = Vl=te-^^ b t- b ^-t Icn 
V = ^2(\Vk + f k )(b{ + b k ), 

k<ui c 



(5) 
(6) 



have zero thermal expectation values, (Vi)T = (i = 
±, 0). Following Silbey and Harris [14], we compute the 
Bogoliubov-Feynman upper bound of the free energy 
[18]: 

-0ii o 



logTrte-"" ) + {Hresh + 0({7&.)t) 



A b = -I: 

(7) 

A first relevant difference with the spin-boson model, is 
that now the Hamiltonian Ho with which the thermal 
averages are taken, contains a potential term quadratic 
in the position of the particle, with potential strength 
g > 0. The effect of such a term is a certain localiza- 
tion of the particle at the center of the chain, making 
the Hamiltonian not easy to diagonalize, as the actual 
eigenenergies and eigenstates depends on t lan and g. 
While a closed expression for Ab cannot be worked out, 
the condition for a minimum does give rise to the fol- 
lowing self-consistent equation for the variational pa- 
rameter f k : 



fk = 



k 



i%^coth(f) 



(8) 



where 7i ,kin = -*ren J2m(\ m ) ( m + 1 1 + H.c). Plug- 
ging this expression into the definitions of the renor- 
malized hopping and the potential strength, we can 
get self-consistent equations for the actual parameters 
of the effective Hamiltonian. If we define ST(t T en,g) = 



1 (Ho,ki, 

2~ 



B(t T . 



, we obtain 



s(t rcn ) 



'/' 

Jo 



duj- 



du; 



tren — tC 
CJCOth(^) 

{lo-st coth(^)] 2 



s T coth(^) 

2 I 



St coth( 



(9) 
(10) 



The study of the variational solution requires the nu- 
merical solution of these two equations. Starting from 
a seed, Hamiltonian Ho is diagonalized to calculate 
the thermal expectation values and the corresponding 
renormalized hopping and potential strength. 
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As pointed out in previous works, the variational 
solution regulates the strength of the coupling between 
bath modes and the particle in a way that preserves 
the spirit of the adiabatic approach [10], but at the 
same time it takes control of the infrared divergences 
introduced by the slow modes of the bath. This sort of 
separation between fast and slow modes is one of the 
main virtues of the approach. 

3.1. Zero temperature 

A further insight into the solution can be obtained 
at zero temperature, i.e., when /3 — > oo. Then the equa- 
tions for the Hamiltonian can be integrated analyti- 
cally, yielding: 



t 1 + 



9 = 



e "0+ 1 

2aso 
so + 1 



(11) 
(12) 



where we have taken lu c = 1 for simplicity. The free 
energy bound reduces to an energy bound, which can 
also be worked out analytically: 



A B (T = Q)=E b = <ft ,kin)o + g(q 2 ) 



(13) 



For every a and t we have a trivial solution of the 
equations, corresponding to t Tcn = and, by virtue of 
their relation, g = 0. This would correspond to a parti- 
cle localized at any site of the chain with a cloud of os- 
cillators dressing it. An interesting implication comes 
from the relation between g and t rcn . As g controls the 
non-homogeneous effects induced by the environment 
in the density profile, no localized solution coming from 
the variational approach will show preference for a par- 
ticular site in the chain, as we mentioned it happened 
in the numerical solution. We will discuss this later. 

The energy associated to the trivial solution is zero. 
The existence of non-trivial solutions of the equations 
with lower energy determines the phase diagram at 
T — 0, in terms of the number of sites M and the pa- 
rameters t and a. As shown in Fig. 1, the region where 
there are non-trivial solutions with lower energy de- 
pends on the number of sites being even or odd. For 
odd chains there is a non-trivial solution with finite i ren 
and g for every value of a. In this solution, the poten- 
tial induced by the environment confines increasingly 
the particle at the center, having in fact a crossover to 
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Fig. 1. Graphical solution of the self-consistent 
equation (11). The function plotted is defined 
as <p(t Tea ,a) = t rcn (l + l/s ) Q e- Q /( s o+ 1 ), where 
s — «o (*ron). The condition for non-trivial solutions is 
thus (f>(t Tcn ,a) = t. In the figures we have chosen two 
representative cases of odd and even chains, M = 5, 6. The 
value of the hopping is t = 0.01. 

a sector in which the particle at all instances is local- 
ized at the center (Fig. 2). Such a crossover is also re- 
flected in the renormalized hopping (inset of Fig. 2) , as 
in this region the hopping gets less renormalized with 
increasing coupling to the bath. For even chains there 
is no longer a non-trivial solution for a > 1, and a 
phase transition occurs to a localized phase in which 
the parity is broken. In the delocalized phase (a < 1), 
the presence of the environment results also in a cer- 
tain confinement of the particle to the two central sites 
in this case. The hopping is further renormalized as the 
coupling increases, an no crossover is observed. We will 
discuss this effect in the last section of the paper. 

3.2. Finite temperature 

For finite temperatures, in general, the self- 
consistent equations have to be integrated and solved 
numerically. Again, the results depend on the number 
of sites of the chain being even or odd. In both cases, 
for couplings a < 1, we find a phase transition at some 
critical temperature T* from a coherent regime where 
trcn and g are finite, to an incoherent high-temperature 
region where t Icn = g = 0. In this phase, the Hamilto- 
nian has degenerate energy levels and they are equally 
occupied in a statistical thermal mixture. In Fig. 3 we 
see a typical plot of the behavior of t Te n{T) as a func- 
tion of the temperature. Notice that for small coupling 
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Fig. 2. Mean squared position deviation for the dissipative 
confined particle in chains of M = 5 and 6 sites. In both 
cases the coupling to the environment leads to a narrowing 
of the density distribution of the particle towards the cen- 
ter. Here, {g 2 )o = 0.25 is the mean squared position corre- 
sponding to a two-site system. Notice that, for odd chains, 
there is a crossover to a region in which the particle is ef- 
fectively localized at the central site. Inset: renormalization 
of the hopping, t Tcn , as a function of the coupling to the 
environment, also for M = 5 and M = 6. The crossover for 
odd chains is reflected in the behaviour of the renormalized 
hopping. 

to the bath, a«l, and high temperatures, T ;§> t Icn , 
we have a trivial expression for the parameter 



so 



1 {T~io.kin)o 

2 {q 2 )o 



M/3t 



2 

rcn 



EM, 
[m - 



m 



(14) 



This is the actual dependence found in [17] for the spin- 
boson model in the same regime, generalized trivially 
to the case of M sites by adding the factor M / Y2m ( m— 
mo) 2 . As stated already in this reference, the critical 
temperature can be worked out analytically and has 
the form: 

1/2 

(15) 



T* 



M 



EM , 



m ) 



Another interesting issue is that, at low temperatures, 
the renormalized hopping i ren gets larger, instead of 
smaller (see Fig. 3). This result was also observed in the 
Ohmic and sub-Ohmic spin-boson model [1,19,17], and 
it is associated to the fact that, at low temperatures, 
only slow bosonic modes are excited, which increases 
the hopping (as opposed to the effect of fast modes). 

In the case of even chains, for a > 1 there is no such 
a phase transition, and only the incoherent phase pre- 
vails. For an odd number of sites, the phase transition 
still occurs for arbitrary a, a reminiscence of the lack 
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Fig. 3. Dependence of the renormalized hopping t 

ren On 

the temperature, for two representative cases of odd and 
even chains, M = 5 and M = 6, respectively. For a < 1 in 
both cases there is a critical temperature that separates an 
incoherent high temperature regime, where f. ron = 0, and a 
coherent low temperature one, where f ron is finite. Notice 
that, for low temperatures, the renormalized value of the 
hopping grows with temperature. 
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Fig. 4. Finite temperature phase diagram of the model. 
The behaviour is different for odd and even chains. In both 
cases there are two phases, one where t Tcn = (incoherent) 
and other where t ren is finite (coherent). For small coupling 
a to the bath, the behaviour of the critical temperature 
is independent of the chain being odd or even. For large 
couplings the behaviour is dramatically different, and the 
coherent phase only prevails for a > 1 for chains with an 
odd number of sites. 



of phase transition at T — 0, that we discussed in the 
previous section. A phase diagram with these features 
is shown in Fig. 4 for the two representative cases of 
M = 5, 6. 
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4. Discussion 

The variational calculation predicts a non-trivial 
phase diagram both at zero temperature and for finite 
temperatures. Being a sort of sophisticated mean-field 
calculation, it gives an effective Hamiltonian for a 
quantum particle with renormalized parameters t Tcn 
and g due to the coupling to the bath. Hence, two main 
effects arise from the environment: the renormaliza- 
tion of the hopping, and the confinement at the center 
of the array due to a quadratic potential. The results 
have a dependence on the number of sites in the array, 
with qualitative differences when odd or even chains 
are considered. 

Let us understand those results. At zero tempera- 
ture, there is no phase transition when the number of 
sites is odd. When the number of sites is even, the re- 
sults are remarkably close to those found for the spin- 
boson model. This is not casual. As shown in Fig. 2, due 
to the quadratic potential, the particle is nearly almost 
localized at the two central sites for coupling strengths 
close to a = f. Thus, we have an effective spin-boson 
model, and the phase transition is essentially the same 
found in this model. Such an argument also applies to 
the finite temperature case, where again we find results 
close to those predicted by the variational calculation 
in the spin-boson model. 

In the case of an odd number of sites in the chain the 
situation is very different: here, at zero temperature, 
the quadratic potential confines the particle around the 
single central site for enough large couplings, and no ef- 
fective spin-boson model physics arises. Once the par- 
ticle is effectively located at the center, the renormal- 
ization of the hopping reverse its tendency, as shown in 
Fig. 2, increasing again. An explanation for this effect 
would come from the fact that the bath is not coupled 
to the central site, as here q\mo) = 0. Thus the renor- 
malization of the particle parameters would be smaller 
as the particle is more localized at the center. This be- 
haviour is reflected in the finite temperature phase dia- 
gram: the coherent-incoherent transition being related 
to the renormalization of the coupling by the bath, the 
critical temperature starts to increase again. For very 
large coupling to the environment, it can be seen that 
the renormalized hopping tends to its bare value. 

The results can be compared with the NRG results 
obtained in [9]. Both approaches qualitatively agree in 



the predictions for the mean squared position of the 
particle in the delocalized phase, in the sense that the 
particle becomes localized at the center, though the 
variational calculation shows a faster degree of local- 
ization as a function of the coupling strength than the 
numerical calculation [If]. 

More relevant is the phase transition found, with 
NRG, at a ~ 1, independently of the length of the 
chain, and where the particle becomes localized in pre- 
ferred sites. A result that disagrees with the phase tran- 
sition predicted in the variational calculation. A way to 
reconcile both approaches is to suppose that the effec- 
tive Hamiltonian for the quantum particle is missing 
higher order contributions to the potential generated 
by the environment. Particularly, the simple potential 
= 9iQ 2 + 92Q 4 cou ld explain the numerical results 
just having g\ oc a c — a and g-z > 0. With this depen- 
dence, the potential would have a minimum at the cen- 
ter for a < a c , and two minima symmetrically located 
around the latter when a > a c . Then the nature of 
the phase transition would be more complicated than 
that found in the spin-boson model, as it would require 
the assistance of the potential generated by the cou- 
pling. Besides, this would predict a non-homogeneous 
density pattern in the localized phase, as contrary to 
what is found in the variational calculation. So far, 
however, more complicated generalizations of the vari- 
ational ansatz have failed to provide such an effective 
potential, and further research is needed in order to 
clarify this point. 
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